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Abstract. We present near infrared velocity curves for 12 long period variables (LPVs) in the globular cluster 47 
Tuc (NGC 104). New light curves are also presented for these variables. Results are compared with the period- 
luminosity sequences occupied by the LPVs in the LMC. Sequence C variables (fundamental mode pulsators) 
have larger velocity amplitudes than sequence B variables (first overtone pulsators). We show that, at similar 
luminosities, higher mass loss rates are associated with larger pulsation amplitudes. One variable (VIS) does not 
fit the normal period luminosity sequences and it has an unusually large amount of circumstellar dust, suggesting 
that it has recently undergone a thermal pulse on the AGB. Finally, we report the discovery of three new long 
period variable stars in the cluster core, all previously found to have a large infrared excess. 
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1. Introduction 

In recent years, there have been significant advances in our 
understanding of pulsation in long period variable stars 
(LPVs). Particularly important is the discovery of multi- 
ple period-luminosity relations for AGB stars in the LMC 
and the interpretation of all but one of these relations in 
terms of radial pulsation in low order modes (Wood ct 
al. ll999(l . Another advance is the realization of the impor- 
tance of pulsation in the production of mass loss from red 
giants (e.g. Wood lTOT^ Bowen & Willson lTMTl Hofner et 
al. lTOMIl . 

One problem with the observational study of pulsa- 
tion and mass loss in red giants is that nearly all the stars 
studied so far are field stars. This means that we have no 
way of directly determining the initial or current stellar 
mass. Any relations that are found for the field stars will 
be broadened by the range in mass (age) and metallicity 
existing among the stars. For Galactic stars, but not LMC 
stars, we have the additional problem that the luminosity 
(distance) is generally poorly known. As a result of these 
uncertainties, precise comparison of observed periods, lu- 
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minosities, amplitudes and mass loss rates with theoretical 
calculations is complicated. A way around these problems 
is to observe pulsating red giant stars in star clusters where 
the giants have a common (initial) mass, composition and 
distance. We can then readily compare theoretical and 
observed period-luminosity relations, and we can sec how 
the mass loss rate depends on pulsation properties and 
quantities such as mass, luminosity and metallicity. 

The driving of mass loss by pulsation is a complicated 
process. Levitation of the stellar atmosphere caused by the 
pulsation leads to the formation above the photosphere of 
a cool and dense environment where dust grains form and 
grow efficiently. Radiation pressure on dust grains then 
combines with momentum in shock waves to drive the 
mass loss (e.g. Wood lWT^ Bowen & Willson [T99II Hofner 
et al. ll996|l . It is found that larger pulsation velocity am- 
plitudes and higher luminosities drive a higher mass loss 
rate. One of the main purposes of this study is to observe 
the pulsation velocity amplitudes in a good sample of red 
variables in a star cluster so that a comparison of mass 
loss rate with pulsation velocity amplitude (and luminos- 
ity) can be made. 

The pulsation velocities of LPVs are best measured in 
the near infrared since optical absorption lines measure 
velocities in the outer layers only. These optical velocities 
are nearly always directed inward relative to the stellar 
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center-of-mass (e.g. Wood ni)79f) . The near-infrared lines 
reveal the deeper, larger amplitude pulsation. Monitoring 
of the velocity variations in field variables using infrared 
(1.6 /Ltm) CO lines has been carried out by Hinkle (|1978|l . 
Hinkle et al. H1982II1997|I and Lebzelter et al. if^nnUjl . These 
observations revealed velocity curves with amplitudes be- 
tween 3 and 30km/s. The velocity curves of Miras are 
sawtooth-shaped and of large amplitude, with phases of 
line doubling around maximum light. The typical velocity 
curves of semiregular variables (SRVs) are of smaller am- 
plitude than in Miras and they are rather sinusoidal (i.e. 
continuous with no line doubling), although they can re- 
flect the partly irregular behaviour seen in the light curves. 
A summary is given in Lebzelter & Hinkle H2U(J2|I . 

The globular cluster 47 Tuc (NGC 104) was chosen for 
this study since it contains the richest known collection 
of long period variables (LPVs). Properties of this cluster 
from the literature are summarized in Sect .[21 Along its 
giant branch, four Miras and more than 10 semiregular 
and irregular variables have been detected (Sawyer-Hogg 
119731 Llovd-Evans ll97l|l . Several infrared searches for in- 
dications of mass loss from giant branch stars have also 
been reported (see Sect.OJ. 



2. 47 Tuc and its variables 

47 Tuc is a prototype "metal-rich" globular cluster as 
well as being one of the closest. Recently, Gratton et al. 
H2003|l presented a very detailed study on the fundamen- 
tal parameters of 47 Tuc. They derive a metallicity of 
[Fe/H]= -0.66±0.04 on the Carretta & Gratton lfTM7|l 
scale, which is in agreement with the findings of several 
other groups (see references in Gratton et al.). A metal- 
licity based on the Zinn & West H1984|l scale is given by 
e.g. Briley et al. fT(^ with [Fe/H]=-0.76. 

An accurate distance to the cluster is still a matter 
of some debate. Values in the literature for the distance 
modulus of 47 Tuc scatter by about ±0.2 mag (for an 
overview see Gratton et al. l2003|l . Gratton et al. H2003|l 
find a distance modulus of (m— M)y=13.50±0.08 from 
main sequence fitting. All studies agree that there is very 
little reddening towards 47 Tuc {£{3 - V) =0.024±0.05). 
Based on their estimates for the distance and metallic- 
ity Gratton et al. H2003|l derived an age of 47 Tuc of 
11.2±l.lGyr. Hesser et al. ()1987|l derived a turnoff mass 
of 0.9 Mq, but used an age estimate of 13.5 Gyr. However, 
with the age from Gratton ct al. and the more recent theo- 
retical isochrones from Bertelli et al. p994|l we again come 
to a turnoff mass between 0.86 and 0.9 Mq. 

A large number of photometric measurements in 
the blue, visual and near- and mid-infrared (Lee 119771 
Frogel ingSI Montegriffo et al. lT^ Origha et al. lTMTl 
Ramdani & Jorissen l2001(l establish a very accurate color- 
magnitude diagram of this cluster showing a well defined 
AGB. Luminosities and surface temperatures have been 
derived for the AGB stars (e.g. Whitelock fTOSSjl . 



Light curves of the LPVs VI to V8^ were measured 
by Arp et al. H1963|l . Several additional variables were de- 
tected by Lloyd-Evans & Mcnzies H1973|l . A further study, 
including the variables V3 to V7, V11(=W12), V13 and 
V18, was presented by Fox H1982fl . The light variability of 
these stars will be further discussed in Sect. l4.ll 

3. Measurements of mass loss in 47 Tuc 

There have been many attempts to measure the mass loss 
from AGB stars in globular clusters. Mass loss at the 
end of the red giant branch (RGB, also called the first 
giant branch or FGB) phase has already been proposed 
to explain the observed gaps in the horizontal branch of 
globular clusters (Soker et al. l2001l Schroder & Sedlmayr 
I2001f) . Cohen H1976|l proposed the existence of circumstel- 
lar shells (and hence mass loss) around first giant branch 
stars in globular clusters to explain emission components 
in the Ha line. Bates et al. (|1990f) summarized different 
studies on Ha lines of globular cluster stars and found a 
further indication for circumstellar shells from the NaD 
line profiles. Lyons et al. H1996|l discussed mass motions 
in the atmospheres of 63 red giants from five different 
globular clusters with the help of Ha and NaD lines. 
Core shifts of these lines indicate mass fiow. The lower 
luminosity limit for outfiow from Ha and NaD lines is 
log(L/LQ)=2.5 and 2.9, respectively. Only part of the red 
giants showed an outflow. Standard deviations in the final 
masses of white dwarfs in globular clusters of the order of 
0.1 mag may indicate a stochastic nature of the mass loss 
(e.g. see the discussion by Alves et al. l20U0|l . 

Not only is mass loss known to be essential in stel- 
lar evolution but the mass loss from red giants in glob- 
ular clusters should result in intracluster gas and dust 
(Evans et al. l2003|l . Searches for the intracluster gas have 
been limited in their results. The low metallicity no doubt 
plays a role in reducing silicate emission, making dust de- 
tection difficult (Frogel & Efias lTM^ Helling et al. lM)^ . 
However, findings from AGB stars in the Large Magellanic 
Cloud (LMC) confirm that low metallicity stars reach 
high mass loss rates (e.g. Wood et al. ll992l van Loon et 
aLMI). 

A few recent papers provide more detailed information 
on circumstellar shells around AGB stars in 47 Tuc. One 
of the most definitive results is by Ramdani & Jorissen 
(|2001|l who used the ISO satellite to measure the mid- 
infrared emission around 6 AGB variables in 47 Tuc. Their 
findings, however, highlight the difficulty of correlating 
pulsational properties with mass loss. Half of the vari- 
ables, namely V5, V7 and V15, show no or only marginal 
12/im-excess while the other three stars, V3, VI 1 and V18, 
do have a detectable excess flux at 12/im. The excess of 
VI 1 is rather small. V3 is a Mira and would be expected to 
show an excess. But the largest infrared excess was mea- 
sured for V18, an irregular AGB variable of only modest 

^ Throughout this paper we use the nomenclature from the 
catalogue of globular cluster variables by Clement et al. II2U0H . 
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luminosity. Origlia et al. (|2UU2|I suggested that the mass 
loss from red giant stars is episodic. They found an IR 
excess for V8 and no excess for V21 (see also Sect. l5.3() . 

Overall, some indication of the presence or absence of 
circumstellar material exists for 12 AGB variables in 47 
Tuc. A substantial infrared excess in V3, Vll and V18 was 
detected by Gillett et al. H1988|l . but these authors did not 
find indications for an infrared excess in V5 and VI 3. An 
excess V-band polarization was found in VI to V4, V6, 
V8 and Vll by Forte et al. H20()2|l indicating the exis- 
tence of circumstellar material around these stars. Glass 
Sz Feast (|1973f) reported an L— band excess in V3. Frogel 
& Elias H1988|l detected a lOfim excess in the four clus- 
ter variables VI, V2, V3 and V4. They calculated total 
mass loss rates for these objects ranging between 4.6 and 
12.3x10"^ MQyr~^. A 10/^m excess was also detected in 
VI and V3 by Origlia et al. (|1997fl . However their esti- 
mated dust mass loss rates are about ten times less than 
those of Frogel & EUas ^E^. 

4. Observations and results 

We selected 12 long period variables in the globular clus- 
ter 47 Tuc for which a period or at least some informa- 
tion on the variability type was given in the catalogue of 
variable stars in globular clusters published by Clement 
et al. I)2(JU1|I . These stars are moderately bright infrared 
sources, with K^6-8. Table ^ summarizes the properties 
of the stars. 

Time series of infrared spectra in the H band were 
obtained in 2001 and 2002 with the 74inch telescope at 
Mount Stromlo Observatory, AustraUa. The NICMASS 
detector, already successfully used for a preceding pro- 
gram at Kitt Peak (Joyce et al. 11998(1 . was used at the 
Coude focus of the telescope. The standard infrared ob- 
servation technique was used. Spectra of each star were 
obtained at two different slit positions to allow sky sub- 
traction. Resolution was set to R=37000. We achieved 
a S/N ratio of 30 or better. An example spectrum is 
shown in the left part of Fig.Q] with an identification of 
some relevant spectral features. The spectral range cov- 
ered a number of second overtone CO lines, some OH 
lines, and a few metallic lines. When 47 Tuc was visible 
at Mount Stromlo spectra were obtained approximately 
once a month. The observing program had an unexpected 
end when Mount Stromlo Observatory was destroyed by 
a bush fire in early 2003. Nevertheless time coverage and 
sampling of the 47 Tuc AGB variables, combined with a 
few additional spectra taken in the same wavelength re- 
gion with the PHOENIX spectrograph at Gemini South, 
is sufficient for the investigation presented here. 

The bright stars a Get and S Oph have been used as 
primary velocity standards (Udry et al. ll999|l . Velocities of 
the variables were determined by a cross correlation tech- 
nique, using the IRAF task fxcor. Typical velocity uncer- 
tainties, determined from multiple observations of some 
stars in the same or consecutive nights, were found to be 
~0.4 kms-^ 
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Fig. 1. Left panel: Example spectrum of the star Vll. 
Several prominent features are identified. Right panel: 
Example of hne doubling (CO 4-1, V3). 



To determine proper phases for the velocity curve, 
photometric measurements of all variables in our sam- 
ple except V5 (see below) were done as well. Data in 
the blue and red MACHO filters were obtained with the 
50inch telescope at Mount Stromlo (also destroyed in the 
January 2003 fires). V and / measurements were ob- 
tained with ANDICAM at the YALO telescope (before 
the MACHO dataset), with a few additional measure- 
ments from ANDICAM at the CTIO 1.3m telescope (after 
the MACHO dataset). The different data sets have been 
combined using transforms based on about 40 nonvari- 
able cluster stars ranging 'm V — I between 0.9 and 2.1. 
The MACHO blue filter has a mean wavelength close to 
that of the V filter so MACHO blue magnitudes can be 
reliably transformed to V. Long time series are therefore 
available in the V filter. No obvious phase shift relative to 
the other filters was noticed, but a conclusive result can- 
not be drawn from our data. From the V light curves we 
determined the time of the light maxima, phase zero. 

4.1. Light curves 

Light curves for most of the variables are shown in the 
upper panels of Figs.|21 to ^1 Periods derived from our 
measurements were in most cases in good agreement with 
the values from the literature. In our data V4 shows a 
main period of '^170 d in agreement with the value of 
165 d from Arp et al. (fBMnjl . Fox iflM^ found a period 
of 82 days for V4 and suggested that this star switched 
between two modes. This type of multimode behaviour 
has subsequently been found to be common in the light 
curves of LPVs (e.g. Wood et al. ll999|l . For V5, we un- 
fortunately do not have sufficient usable data as this star 
was located in a bad area of the CCD. The Sawyer-Hogg 
value H1973|l is 60 days, but different values can be found 
in the literature, ranging between 29 and 70 days (Fox 
I1982|l . From our velocity data we rather favour a period 
of 50 days. The short periods of 52 and 100 days reported 
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Table 1. Properties of LPVs in 47 Tuc. 



Name a (2000) 5 (2000) J K P [d] 



VI 


00 


24 


12.4 


-72 


06 


39 


7.45 


6.21 


221 


V2 


00 


24 


1S.4 


-72 


07 


59 


7.52 


6.29 


203 


V3 


00 


25 


15.9 


-72 


03 


54 


7.49 


6.27 


192 


V4 


00 


24 


00.3 


-72 


07 


26 


7.S7 


6.69 


165,82 


V5 


00 


25 


03.7 


-72 


09 


31 


8.65 


7.47 


50 


V6 


00 


24 


25.5 


-72 


06 


30 


8.54 


7.43 


48 


V7 


00 


25 


20.6 


-72 


06 


40 


8.18 


6.97 


52 


V8 


00 


24 


0S.3 


-72 


03 


54 


7.94 


6.70 


155 


Vll 


00 


25 


09.0 


-72 


02 


17 


7.91 


6.71 


52,160? 


V13 


00 


22 


5S.3 


-72 


06 


56 


8.79 


7.70 


40 


VIS 


00 


25 


09.2 


-72 


02 


39 


8.59 


7.47 


83? 


V21 


00 


23 


50.1 


-72 


05 


50 


8.07 


6.78 


76 



Notes: For V18 and V21 both period determinations are from 
this study. For VI we suggest a period of 221 days from our 
data instead of 212 days listed in the literature. For V5 we 
favour a period of 50 days (instead of 60 days listed in the 
Clement et al. catalogue). For Vll see text. All other periods 
are from the catalogue of Clement et al. I)2()()l|l . The J and 
K magnitudes are averages of the maximum and minimum 
observed values from the following sources: the 2MASS Point 
Source Catalog; Fox H1982|l : Menzies & Whitelock l|1985^ : and 
Frogel et al. I)1988|l . All magnitudes have been converted to 
the AAO system using conversions in Allen & Cragg (I1983II 
and Carpenter Ij2004|l . 

for Vll could not be confirmed in our data. Instead, the 
star shows some long period variation lasting more than 
200 days fFig. llO() . Compared to the data presented by Fox 
l|1982|l the amplitude of the variation is smaller. This star 
appears to be another multimode pulsator. For V13, there 
is a hint of the catalogued ^40 d period in our light curve, 
but a period of ~90d may also be present. Additionally, 
our light curve suggests a long period fFig. lll|l about 10 
times the short period. This star appears to be yet another 
multimode pulsator. 

V18 was reported variable both in the catalogue of 
Sawyer-Hogg H1973|l and by Lloyd Evans (|1974(l , while Fox 
l|1982|l found no variability in this star. We found photo- 
metric variability with an amplitude of about 0.2 mag in 
this object (Fig. 112(1 . Our short light curve would suggest 
a period of about 83 days. For V21, no period determina- 
tion existed in the literature. Our data set (see Fig. I13|l 
shows a period of 76 days. 

4.2. Velocity curves 

Velocity curves were determined for all 12 stars of our 
sample. All stars show at least some velocity variability 
above our detection threshold. Line doubling was detected 
in three stars of our sample (Tab. 01. An example of line 
doubling is shown in the right part of Fig.^ 

The velocity curves we derived are shown in the lower 
parts of Figs.Eltoll3ltogether with the corresponding light 
change. The velocity curves can be roughly separated into 
three groups: VI, V2, V3, V4, and V8 show velocity curves 
very similar in shape to those typical of Miras found in the 



Table 2. Data on the AGB variables in 47 Tuc. 



Name 


velocity amplitude 
[km/s] 


variability 
type 


comments 


VI 


20 


Mira-like 


line doubling 


V2 


23 


Mira-like 


line doubling 


V3 


22 


Mira-like 


line doubling 


V4 


18 


Mira-like 


double-period 


V5 


8 


regular 




V6 


7 


semiregular 




V7 


4 


semiregular 




V8 


16 


Mira-like 




Vll 


4 


irregular 




V13 


12 


long period 


double-period 


VIS 


5 


irregular 




V21 


7 


long period 





solar neighborhood (Lebzelter & Hinkle 2001). The veloc- 
ity amplitude of VI, V2 and V3 are all similar to nearby 
Miras, too. The other two stars show a similar shape but 
a clearly smaller amplitude. While we confirm the Mira- 
like nature of these five variables (e.g. noted by Whitelock 
I1986|) . we note that the smaller velocity amplitude of the 
latter two stars is also accompanied by a smaller light am- 
plitude. The second group of variables, consisting of V5, 
V6 and V7 is comparable to local semiregular variables 
(Lebzelter & Hinkle 2001). Their amplitudes are much 
smaller than in the first group. V5 shows a very regu- 
lar light change, while the other two objects are obviously 
not strictly periodic, in agreement with the semiregular 
nature of their light curves. The third group consists of 
Vll, V13, V18, and V21. In all four of these stars, there 
is no obvious correlation between the light change and the 
velocity variations. In the case of V13, it cannot be ruled 
out that variations occur on a time scale similar to the 
short or even the long period (see above). For Vll our ob- 
servations do not allow a clear picture of the light change, 
so a correlation is not possible. On the other hand, V18 
and V21 do show a rather well defined light change, but 
velocity variability occurs on a different time scale. The 
reason for that is not clear and longer spectroscopic and 
photometric time series will be needed to understand this 
phenomenon. Table [3 summarizes the results on the ve- 
locity variations in the 47 Tuc AGB variables. The table 
gives the total velocity amplitude, the characteristics of 
the velocity curve and comments on the occurrence of line 
doubling. 

5. Discussion 

5.1. Pulsation and Stellar Evolution 

Observations from the Magellanic Clouds indicate that 
AGB variables are found on four distinct period- 
luminosity relations (e.g. Wood l^OOfl . Three of these se- 
quences can be interpreted in terms of stars pulsating in 
the fundamental mode and the first, second and third over- 
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Fig. 2. Light (upper panel) and velocity variations (lower 
panel) for VI. Filled symbols in the lower panel indicate 
individual velocity measurements. For a better illustration 
of the velocity change data are repeated shifted by an 
integral number of periods forward and backward in time 
(open symbols). The period listed in Table 1 is used. The 
typical error bar for the velocity data is indicated. 
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Fig. 4. Same as Fig.^lfor V3. 



tone modes^. The fourth period- luminosity relation still 
lacks a definite interpretation (Wood et al. l2(JU4fl . For the 
Milky Way AGB variables the absence of reliable distances 
limits the establishment of such a period-luminosity rela- 
tion (e.g. Bedding & Ziilstra [l998|l . In the case of both the 
LMC and the Milky Way, the diversity of the stars in age, 
and therefore initial mass, complicates the interpretation. 

The AGBs of globular clusters are characterized by a 
much more homogeneous set of stellar parameters. The 
richness of 47 Tuc in AGB variables allows us to use this 
advantage for a discussion of stellar pulsation in relation 
to luminosity, i.e. evolutionary stage. Period-luminosity 
diagrams for globular clusters have been constructed by 
several authors (e.g. Feast et al. l2l)U2|l . Typically, only a 
relation for the large amplitude, long period variables 
('Miras') was derived. This relation is in agreement with 
the sequence C found from LMC data (Wood 0000). In 
Fig. ^3 we placed the AGB variables of 47 Tuc studied in 
this paper into a logP-K-diagram. The approximate loca- 
tion of the P-K-relations B (first overtone) and C (funda- 
mental mode) determined in the LMC fWood are 
indicated by lines. Relations were transformed from the 
LMC to 47 Tuc distance using a LMC distance module 
of 18.515 (Clementini et al. l2mJ3|l . Stars with more than 
one well established period in the literature or from our 
own light curve data are plotted twice and are marked 
accordingly in Table |21 



Fig. 3. Same as Fig.Hfor V2. 



^ Sequence A is attributed to both second and third overtone 
mode pulsators. 
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Fig. 7. Light (upper panel) and velocity (lower panel) 
variations for V6. A typical error bar for the velocity data 
is indicated. 
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Fig. 6. Velocity change of V5. Only part of the time of 
monitoring is shown. Symbols as in the lower panel of 
Fig.[21 No parallel light curve data exist for this star. A 
typical error bar for the velocity data is indicated. 



It can be seen that the variables nicely follow these two 
relations. The good fit of the P-L-rclations to the 47 Tuc 
data after scaling from the LMC is indeed remarkable. 
Obviously, the physical mechanism that determines the 
location of at least these two sequences in the logP-K- 
diagram is the same for LMC and 47 Tuc. Furthermore, 
the distance moduli we used in this paper seem to be very 
well determined at least concerning the relative distance 
between LMC and 47 Tuc. However, it has to be noted that 
the location of the sequences in the LMC is not accurately 



determined, as there is a significant scatter (probably due 
to the larger scatter in stellar parameters and maybe also 
the extension of the LMC) around the average locations 
used here. We therefore do not want to overinterpret our 
findings on that point. 

There is a clear bipartition in our logP-K-diagram 
(Fig. injl : At lower luminosities only sequence B is pop- 
ulated (see below for the case of V18), no variables are 
known in 47 Tuc that would fall on the lower part of logP- 
K-rclation C. It is very likely that the stars V5, V6, V7, 
and VI 3, which we find in this lower luminosity part are 
probably in an earlier evolutionary state than the stars 
at higher luminosity (see below). In a first approach we 
assume that the luminosity of the stars in Fig.^] corre- 
sponds to the evolutionary status. 

At the highest luminosities only sequence C is occu- 
pied, and the stars there are the Miras VI, V2, and V3. 
Below these stars on sequence C lie V4 and V8. Both stars 
are located close to the tip of the RGB. They have similar 
velocity amplitude and K brightness. V4 appears to switch 
modes between sequence C (P = 165 d) and sequence B 
(P = 82 d). From our data we have to conclude that V4 is 
currently only on sequence C. VI 1 having approximately 
the same K luminosity as V4 is another star suspected of 
showing mode switching (see Sect. 14.1(1 located clearly off 
sequence B towards shorter periods (i.e. near sequence A). 
The period used in the plot is taken from Fox H1982|l . It 
is interesting that both stars with possible mode switch- 
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ing arc found at similar K luminosity, namely close to the 
RGB tip. 

For V13 we found indications of a second period, and 
the star is therefore shown twice in Fig. El The long pe- 
riod of V13 is located close to sequence D from Wood 
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()2()()()(l . The nature of this period- luminosity relation is 
not yet understood (Wood et al. l2UU4|l . 

Symbol size in our P-K-diagram is linearly propor- 
tional to the measured velocity amplitude. In the case of 
V4 we stress again that the shorter period, marked with 
an open symbol in Fig. 1141 is clearly not dominating the 
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light or velocity variations at the moment. The symbol 
size of the V4 data point on sequence B may therefore 
be misleading as we have no indications that the veloc- 
ity amplitude would be the same if the 82 days period 
would dominate. According to our data the star is cur- 



log P 

Fig. 14. logP vs. K diagram for the long period vari- 
ables in 47 Tuc. Symbol size denotes the velocity ampli- 
tude (ranging between 4 and 23kms~^). Lines indicate 
the approximate location of sequences B and C found by 
Wood (|2000|l for LMC long period variables shifted to the 
distance of 47 Tuc. The dotted horizontal line marks the 
tip of the RGB according to Ferraro et al. H2()()()|l . Stars 
with possible multiple periods arc shown twice. 



rcntly monoperiodic. It can be seen that the large ampli- 
tude variables arc all found along sequence C. Our results 
indicate further that the velocity amplitude is increasing 
along sequence C. Velocity amplitudes along sequence B 
are much smaller and do not show a steady increase of ve- 
locity amplitude with luminosity. The large amplitude of 
V13 seems to be associated with the long period (^400 d) 
variation rather than the 40 d variation of sequence B. 

These results lead us to the following probable evolu- 
tionary scenario: Stars evolve up the AGB from low lumi- 
nosities. They pulsate first in the first overtone (sequence 
B) then, at intermediate luminosities close to the RGB tip, 
they go through an interval of mode switching back and 
forth between first overtone and fundamental mode and 
finally, at the highest luminosities, they remain pulsating 
in the fundamental mode. There is clearly an increase in 
both the light and velocity amplitude associated with the 
switch from first overtone to fundamental. The large veloc- 
ity amplitudes found at the tip of sequence C are similar 
to what is expected from models for fundamental mode 
pulsation (e.g. Scholz & Wood 12101111 • 
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The role of V18 in this scenario is not clear. As noted in 
Section|21 this star has one of the largest infrared excesses 
among the 47 Tuc variables. In Fig.^lwc find the star in 
between sequences B and C at a rather low luminosity. Its 
velocity amplitude is rather small. However, locating this 
star on the period axis is difficult as the star has shown 
periodic, nonperiodic and constant phases in the past. 

Comparing our results with evolutionary models re- 
quires some adaption of this scenario: Due to luminosity 
variations during a Thermal Pulse the current location of a 
star in the logP-K-diagram may not correspond to its evo- 
lutionary status along the AGB. However, as was shown 
e.g. by Boothroyd & Sackmann H1988|l the probability to 
observe the star at luminosity values highly deviating from 
the mean (and for the star's evolutionary status typical) 
value is comparatively low. This fact is even more ex- 
pressed in stars of 1 Mq or less where the intcrpulsc time 
is very large compared to the duration of the pulse itself 
(Vassiliadis & Wood fn)93l) . From the models of Vassiliadis 
& Wood (|T^ see their Fig. 13) for a star of IMq and 
Z=0.004 we would expect to see about 1/6 of our sample 
in a phase strongly deviating from the mean global pa- 
rameters of the star, in our case these would be 2 stars 
(if our sample consists of AGB stars only, see Sect. I5.2|l . 
Most probably, these stars are expected to be found at a 
lower luminosity and a shorter period. It is therefore un- 
likely that one of the stars on sequence C belongs to this 
group. Following these considerations we have to add to 
our scenario the possible occurrence of loops in the logP- 
K-diagram during the evolution up the AGB. VIS with its 
large IR excess and odd position in the logP-K-diagram 
may be an example of such a star. However, it is unlikely 
that e.g. the mode switch of V4 is an indication of such a 
loop because of the time scale of the mode switch. Arp et 
al. H1963(l gave a period of 165 days from data obtained 
in 1955/1956. Our light curve also favours this period. In 
between these two measurements are the light curve data 
presented by Fox H1982ii giving the shorter period. The 
star therefore switched its pulsation mode back and forth 
within about 50 years. This seems to be hardly compati- 
ble with the longer timescales expected from evolutionary 
models. The mode switches seem to be concentrated in 
the luminosity range close to the RGB tip. 

5.2. AGB or RGB stars? 

Throughout this paper we have assumed that all the vari- 
ables discussed here are on the AGB. As can be seen from 
Fig-El a significant part of our sample is located below 
the tip of the red giant branch (RGB). The RGB tip plot- 
ted in Fig. ll4l is taken from the recent study of Fcrraro 
et al. pHinjl . They give a value of KTi?,GB=6.75±0.2mag. 
This result is in good agreement with the theoretical value 
expected for the mctallicity of 47 Tuc (Salaris & Cassisi 
I1998|l . All stars are above the low luminosity limit for 
thermally pulsing stars, so that they could belong to ei- 
ther the RGB or the AGB. The stars well above the RGB 



tip are clearly AGB stars, namely VI, V2, V3. V4 and 
V8 are both fundamental mode pulsators slightly above 
the RGB tip and thus they can be counted as AGB stars, 
too. For the other stars of our sample, we cannot decide 
on their AGB nature from luminosity alone. Colour infor- 
mation may be a good indicator for a separation, but our 
photometry is not homogeneous enough to decide on this 
rather small effect as a mean colour over the light cycle 
would be required. 

Can variability properties allow us to distinguish be- 
tween AGB and RGB stars? Using artificial luminosity 
functions, Alves et al. p998|l and Wood et al. (|1999|l ar- 
gued that most pulsating stars below the RGB tip are 
in fact AGB stars. However, using data more sensitive to 
smaller amplitude variables, Ita et al. (|2()()2(l and Kiss & 
Bedding l|2()()3|l have shown that a substantial fraction of 
variables below the RGB tip should be RGB stars. 

According to the logP-K-diagram for LMC variables 
given by Kiss & Bedding <|2()()HI see their Fig. 4), most of 
the variables detected below the tip of the RGB are pul- 
sating in the second or third overtone mode. Furthermore, 
Kiss & Bedding split their logP-K-diagram into six parts 
according to the light amplitude in /. Most of the vari- 
ables below the tip of the RGB show amplitudes between 
0.01 (their detection limit) and 0.14 mag. The luminos- 
ity function for these stars indicates a substantial frac- 
tion arc on the RGB. However, for larger light ampli- 
tudes, the RGB component is no longer detectable. Kiss & 
Bedding also give colour information showing that by far 
the largest fraction of stars below the RGB tip with ampli- 
tudes >0.14mag is found on the blue side of the red giant 
branch, which again favours attributing these stars to the 
AGB. Thus larger light amplitudes seem an indication of 
AGB status. 

Kiss & Bedding l|20()3|l used / band data, whereas we 
have V band light curves. Fox H1982|l presents V and / 
light curves for several of our sample stars. For the stars 
of lower luminosity like V5 or V6, the / amplitude is about 
a factor of 2 less than the V amplitude. The V band am- 
plitudes of the stars of our sample arc all >0.2 mag, i.e. we 
estimate an / amplitude of > 0.1 mag. 

V21, V5 and V6 all show light amplitudes of several 
tenth of a magnitude. Such amplitude values favour an 
AGB nature for these stars as well. V7 shows a rather 
small light amplitude, and for VI 1 the parameters of the 
light change are not clear. Therefore, we cannot decide if 
these two objects are RGB or AGB stars. Interestingly, 
both objects also have smaller velocity amplitudes. 

Another unclear case is V13. The total amplitude mea- 
sured by us is Ri0.7mag. But this amplitude is dominated 
by the long period variation, the short period variation has 
a much smaller amplitude. The amplitudes used in the dis- 
cussion by Kiss & Bedding are derived from Fourier analy- 
sis of the individual frequency components. Therefore the 
amplitude of the short period is the relevant one in our 
argumentation. As a result V13 may well be a RGB star. 

Our approach is based on statistical arguments and 
not on a complete understanding of RGB pulsation. In 
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particular, we note that the LMC stars of Kiss & Bedding 
(|2003|l are likely to be more massive than the 47 Tuc stars 
and their amplitude behaviour may be quite different to 
that of 47 Tuc stars. We therefore cannot exclude the pos- 
sibility that some of the variables below the RGB tip are 
on the RGB, but for many of them a location on the AGB 
seems to be more likely. 

5.3. Comparison with mass loss 

Using the velocity amplitude we can separate the logP- 
K-diagram (Fig. ll4|l of 47 Tuc into three to four regions. 
At highest luminosity we find the Mira variables (VI, V2 
and V3) with the largest velocity amplitudes. Going down 
sequence C, wc find the second group (V4 and V8) with 
intermediate amplitudes. The third group consists of stars 
with low velocity amplitude on sequence B (and possibly 
A). Finally, V13 may form a special case with its large 
amplitude due to a long secondary period. 

Mass loss data available in the literature are limited 
to measurements of circumstellar dust. As there may be 
a dust-free mass loss as well, e.g. during the RGB phase, 
we cannot rule out that stars with indications for no or 
only small amounts of circumstellar dust still have a con- 
siderable mass loss rate. However, it is very likely that the 
existence of dust enhances the mass loss rate (e.g. Hofner 
et al. ll996|l . so we expect that the infrared excess is at least 
a measure for the relative mass loss rate within the cluster. 
As a caveat we have to mention the lack of infrared pho- 
tometry beyond 12 /um for the stars of our sample. This 
certainly limits the reliability of estimates for the amount 
of dust surrounding a star. 

Note that the circumstellar material currently around 
each star may originate from a time when the star's pulsa- 
tion was different. Another factor to keep in mind is that 
it is well known (e.g. Olofsson et al. l2000|l that, at least in 
some cases, mass loss in AGB stars may be episodic rather 
than continuous. This has also been noted in the context 
of globular clusters by Origlia et al. (|2QD2J. 

According to Reimers' Hiy75|l law, we would expect the 
mass loss rate to increase with luminosity. Therefore it is 
difficult to separate the effects of luminosity and pulsa- 
tion (velocity amplitude) on the mass loss rate: high mass 
loss rates arc found at high luminosities, where we also 
find large velocity amplitudes. We can, however, try to 
compare mass loss rate and velocity amplitude for some 
stars of similar luminosity. Such a group of stars is found 
around the RGB tip (V8, V4, Vll and V21). Origiia et 
al. measure A'-[12] values for V8 and V21. V8 

shows a clear dust excess while V21 does not. Comparing 
the velocity amplitudes for these two stars shows that V8 
also has a significantly larger velocity amplitude than V21. 
V4 is also a star with an infrared dust excess (Frogel & 
Elias I1988|l . It has a similar velocity amplitude as V8. 
Gillett et al. H1988|l report an infrared excess for Vll, but 
a separation between Vll and the nearby variable V18 is 
very difficult on the IRAS images they used. Ramdani & 



Jorissen ()2(J(J1|I later found an infrared excess of V18 from 
ISO observations, while Vll has only a very small infrared 
excess. We conclude that, at similar luminosity, stars with 
higher pulsation amplitudes show higher mass loss rates. 

V5, V6 and VI 8 form a second group of variables at 
similar luminosity. Unlike the group of stars around the 
RGB tip, these stars all have very similar velocity ampli- 
tudes. Neither V5 nor V6 has an infrared excess reported 
in the literature. But V18 has an outstanding infrared ex- 
cess according to Ramdani & Jorissen ()2001|) . As noted 
above, the high infrared excess of V18 at a relatively low 
luminosity suggests that this star is currently in the lumi- 
nosity minimum following a thermal pulse on the AGB. 

V13, the star with the long secondary period, shows 
no indication of circumstellar dust. It appears that these 
long periods on sequence D are not directly related to a 
mass loss phenomenon (compare also Wood et al. l2UU4|l . 

Origlia et al. H2()()2|l detected five further stars nearer to 
the core of 47 Tuc that show considerable infrared excesses 
but that are not known to be variable star. We checked our 
photometric monitoring data covering two months in late 
2003 for variable stars at the positions given by Origlia 
et al. Four of the five stars had variable counterparts on 
our V frames. The fifth object is the one with the lowest 
mass loss rate in the list of Origlia et al. All four ob- 
jects varied during the two months by several tenths of a 
magnitude. More detailed results on these stars including 
period determination will be presented in a forthcoming 
paper (Lebzelter & Wood in prep.). 

Summarizing our results regarding the correlation of 
mass loss and velocity amplitude: Stars with significant 
mass loss are all pulsating, a fact already noted by other 
authors (e.g. Ramdani & Jorissen l!^001(l . Stars on sequence 
C combine large velocity amplitudes with comparably 
high mass loss rates. On sequence B all stars have low ve- 
locity amplitudes and most of them also have no or very 
modest mass loss. For stars of similar luminosity, increased 
pulsational amplitude seems to significantly increase the 
mass loss rate. One star (V18) currently showing small 
pulsation amplitude has considerable amounts of circum- 
stellar dust, presumably from a past interval of higher 
amplitude pulsation. We speculated that V18 is currently 
doing a loop in the logP-K-diagram as a consequence of 
a recent or ongoing thermal pulse. Its circumstellar dust 
may then originate not only from a phase of higher am- 
plitude pulsation but also from a time the star was more 
luminous. 
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